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The Fourier-Bessel expansion of a function on a circular disc yields a simple series representa- 
tion for the end-to-end probability distribution function w{R, (f)) encountered in a planar persistent 
random walk, where the direction taken in a step depends on the relative orientation towards the 
preceding step. For all but the shortest walks, the proposed method provides a rapidly converging, 
numerically stable algorithm that is particularly useful for the precise study of intermediate-size 
chains that have not yet approached the diffusion limit. As a practical application, we examine the 
force-extension diagram of various planar polymer chains. With increasing joint stiffness, a marked 
transition from rubber-like behaviour to a form of elastic response resembling that of a flexible rod 
is observed. 
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I. INTRODUCTION 



Today, the concept of the random walk is an ubiq- 
uitous tool in the physical sciences, and a multitude of 
different mathematical and physical problems have been 
subsumed under this catch-all phrase JJ |2| ■ Originally, 
the term has been coined by Pearson yf who was inter- 
ested in the probability distribution wp{R) for the end- 
to-end distance in a chain of N randomly chosen steps of 
equal length I in a plane. This article concerns a random 
walk process that is closely related to the initial formu- 
lation of the problem: While Pearson's general setup is 
maintained, we introduce an angular bias into the rela- 
tive orientation of successive steps that is governed by 
a probability distribution p{9). Known as the continu- 
ous persistent random walk in the plane, this modifica- 
tion continues to attract considerable theoretical interest 
1^ ^, i5„ lZJ . Practical applications of the planar walk 
cover a wide range including light scattering, where sums 
of random complex numbers are studi ed M , its related 
use in crystal structure analysis 0, 0, Il0|. and surface 
diffusion processes. The model of the persistent walk ex- 
amined here is of particular interest to the configurational 
statistics of polymers 0| , where the inherent stiffness of 
chemical bonds restricts the relative orientation of adja- 
cent monomers, and obviously for two-dimensional loco- 
motion problems in biology and medicine 0, 0, 0, 0| ■ 

In several of these publications d, IE Q , only mo- 
ments of the probability distribution function w{R, (j)) 
were calculated, while it appeared desirable to the au- 
thors to have access to w{R, (j)) beyond the diffusion limit 
N ^ oo, where this function ultimately takes on a simple 
Gaussian shape. Unfortunately, the search for explicit so- 
lutions even for Pearson's original problem has remained 
largely elusive. While some recent contributors resorted 
to numerical Monte Carlo simulation 0| , in a pioneering 
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work Barakat |l5j | realized that expansion into a discrete 
Fourier series offers a superior alternative in computing 
the radial distribution function w{R) of isotropic random 
walks. The corresponding Fourier-Bessel series solution 
for the Pearson random walk is equally due to Barakat 
0, and subsequently has been advantageously applied 
in crystallographical analysis 0, 0- In this article, we 
extend the Fourier analysis to the persistent random 
walk problem characterized by an anisotropic end-to-end 
distribution function w{R, (j)) and find excellent numeri- 
cal convergence for reasonably smooth angular distribu- 
tion functions p{9), except for the shortest walks. This 
renders the method particularly suited for the study of 
intermediate-size chains whose behaviour conspicuously 
deviates from the diffusion limit. Our derivation of the 
series representation rests on rather elementary methods. 
(A more complicated general analysis based on the char- 
acteristic function c(w, x) has been put forward by Weiss 
and Shmueli When applied to the persistent random 
walk, it implicitly yields the same result.) 



Let us briefly outline the organization of this article: 
In Sectional we develop the theory of the persistent pla- 
nar walk. Starting with moments of the distribution, we 
proceed to derive the characteristic function c{uj, x) m 
form of a matrix product, which we use to deduce an in- 
tegral representation of 'w{R, (j) ) re miniscent of Kluyver's 
solution to Pearson's problem Finally, we establish 
the Fourier-Bessel eigenfunction expansion of this distri- 
bution. Numerical results obtained with our method are 
presented in Section [llll where we study the convergence 
behaviour of the series and display the frequently dis- 
tinctive shapes of the distribution function obtained for 
several choices of the angular weight p{9). Some possi- 
ble extensions of the theory are discussed in Section IIVI 
Finally, a proof of the "projection theorem" that forms 
the basis of our proposed series expansion is sketched in 
Appendix EI while Appendix IbI covers the calculation of 
moments in the Fourier-Bessel scheme. 
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FIG. 1: The geometry of the planar random walk. The planar 
walk is described by the set of relative angles {^i, 62, ■ ■ ■ ,0n} 
formed by adjacent steps, or alternatively by the set of ab- 
solute angles {9i, 62, ■ ■ ■ ,0n} of the segments with respect to 
the a;-axis. (Note that 6\ =61.) 



II. THEORY OF THE 2D PERSISTENT WALK 

We start out with some basic conventions. For sim- 
plicity of notation, in this paper we deal exclusively with 
a planar persistent random walk built of N steps with 
uniform length but the generalization of our theory to 
incorporate individually different step lengths is straight- 
forward. The walk is then fully characterized by the set of 
relative angles {0i,62, ■ ■ ■ , On} at the links of successive 
segments, where the initial angle 0i is measured relative 
to the X-axis of our coordinate system (Figure Q]) . In 
practice, the absolute angles 0^ formed between the v.th 
step and the x-axis as direction of reference are equally 
important, and they form an alternative description of 
the chain. As easily inferred from Figure ^ both sets of 
angles are interrelated via: 



E' 

k=l 



(1) 



Since we are mainly interested in the end-to-end distribu- 
tion function w(i?, (j)), we start the random walk from the 
origin of our coordinate system. We denote its endpoint 
by {Lx,Ly) in cartesian and {R,(f>) in polar coordinates, 
respectively. 

For the persistent walk, the probability of a certain 
configuration p{9i, 62, ■ ■ ■ , On) depends only on the rel- 
ative angles 61, formed by adjacent steps. Assuming 
that these angles are not correlated, and that all links 
are equivalent (including the initial "joint" between the 
first step and the x-axis), this probability reduces to a 
product of independent factors: p{9i, ■ ■ ■ ,0n) = p{Qi) ■ 
... ■ p{9n), where the angular probability distribution 
p{6) = p{9 + 2tt) is a periodic function of the change 
of direction 6. This suggests that a Fourier expansion of 



the function p{9) will be convenient [T^ll^ : 



Pie) 



1 

2^ 
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dOe'^^piO) 



(2) 

Proper normalization of p(0) requires po = 1; the norm of 
all other (generally complex) angular Fourier components 
is then less than unity, \p^\ < 1 for v ^ 0. (Equality may 
occur only for rigidly linked chains with p{9) — S{9 — 60) 
which we exclude here.) With these definitions, we are 
now set to determine quantities of interest for the walk. 



A. Persistence length and diffusion constant 

Our first investigation concerns the average position 
{{Lx)N,{Ly)N) of the chain end and its mean square 
displacement (L^ -I- L^) jv as a function of the number 
of steps N . Both quantities are accessible in a simple 
calculation that introduces some techniques that will be 
commonplace below. It is convenient to interpret vectors 
in the x — y-plane as complex numbers; with this con- 
vention, the vector formed by the :^.th element {l^x \ 4"^' 



) 

is characterized by its absolute angle 9y with respect to 
the a;-axis (Figure : l^x'' + 14"^^ — ' e'^" . The position 
of the chain end is simply the sum of all constituent step 
vectors, Lx + \Ly — I Y^^=i "3'^" i and its properly weighed 
average over all configurations of a A^-element walk cor- 
respondingly reads: 



{Lx+iLy)^=lY, / dM^l) 

N 



AONPieN) e" 



(3) 



l/=l 



Pi 



Here, we successively used the resolution of the absolute 
angle 6*^ into relative angles Ok , the definition of the 
angular Fourier components in |(2Jl, and finally the nor- 
malization condition pq = 1. Wc infer that the average 
displacement of the chain end depends only on the first 
Fourier component of p{0), and if pi vanishes, the distri- 
bution is centered around the origin for all N . Unless the 
chain is absolutely rigid {\pi \ = 1), the mean end-to-end 
vector approaches a well-defined limit as A'^ — > cxd: 



N 



Um {Lx + \Ly) 



VI N 
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Pi 
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Pi 



(4) 



which generalizes the notion of the persistence length of 
the chain put forward in Ref. [TH . For comparison, we 
also state this result using real parameters. Introducing 
Pi = ae'^, where /3 represents the average "slant" angle 
in the angular distribution p{9) one finds in the 
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long-chain limit: 



la{cos j3 ~ a) 
1 ~ 2a cos /3 + 

Zasin/3 
1 — 2a cos f3 + a'^ 



(5) 
(6) 



which yields for the average absolute displacement (R) 
+ (i,)2)i/2 as iV ^ oo: 



la 



\/\ — 2a cos /3 + 



This simple expression might be useful in the interpreta- 
tion of experimental data Q ■ 

The determination of the mean square displacement 
(i?^) TV = (-^x + -^y) N is only slightly more involved. In 
the complex notation, i?^ = [L^ iLy){Lx — ^Ly) = 
l'^ J2^ii=i e''^''~^^'\ and the same series of arguments we 
used in © yields: 



ly.fj,—! 



N 



(8) 



(5=1 



where we exploited 9i, — 0^ = + . . . + 9^ [v > ^) and 
sorted the terms in the double sum according to their 
"distance" 5 — \v — The arithmetic-geometric pro- 
gression in (|SJ) is easily evaluated : 



N 
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(9) 



(A very similar result is obtained for the mean square dis- 
placement in the three-dimensional freely rotating chain 
[nl.) Since the complex conjugate term is formally ob- 
tained by replacing pi with p_i in @, the mean square 
displacement, like the absolute displacement J^l, depends 
only on the first angular Fourier component p±i. For 
Pi = 0, (O reduces to (-R^)^ = NP, the trivial result 
obtained in Pearson's walk. In the limit iV ^ oo, we 
arrive at the following diffusion constant (characteristic 
ratio) Dm = {R'^)n/N for the persistent walk: 



Pil -Pip.i) 



P{l-a^) 



(1 -pi)(l 1 - 2a cos /3 q;2 



(10) 



The asymmetry parameters a and /? of the angular distri- 
bution p{9) thus determine the persistence length {R)oc 
{Tj) and the diffusion constant Doo of the walk and vice 
versa. We point out that for /3 7^ 0, with increasing a 
Doo starts to diminish from a certain stiffness of the joint 
and actually vanishes as a — > 1. This property reflects 
the fact that the walk then tends to move around a cir- 
cle, and no net transport takes place. (An alternative 
derivation of H1U|) has been given by Larralde 0.) We 



conclude this section by noting that higher momenta of 
the distribution function 'w{R, (p) evidently are available 
from similar calculations. However, in practice they are 
more easily evaluated from a Fourier-Bessel series repre- 
sentation (Appendix IB|| . 



B. The characteristic function 



In the previous section, we were able to determine the 



^7) moments {L^ + iLy) m Q and (i?^ 



AT (|8j) of the prob- 
ability distribution w{R, (j)) by summation of some ge- 
ometric series. We now condense the ideas implicitly 
introduced there into a single, coherent formalism and 
employ it towards the corresponding characteristic func- 
tion c(ijj,x) = (exp(ia; • L))jv of the persistent walk. 

So far, we evaluated averages of the form (exp(i0i,)) 
and (exp[i(0i, — ^^)]), respectively. Both expressions may 
be viewed as special cases of the more general average 
(exp(i^^, Ay^i/)), where {Ai, . . . , \n} is a set of N inte- 
gers and the sum covers the range 1 < < A'^. In view 
of the relation 9^ = 9i -\- . . . -\- 9^ ^ between absolute 
and relative angles, its expression in terms of the angular 
Fourier coefficients p\ (0) is straightforward: 



— PA„PAjv+Ajv-i •• • ■•PA„ + ... + A2+Ai ■ 

(11) 

The significance of this development is as follows. Con- 
sider a functional Q{9i, . . . ,9^) of the random walk that 
depends only on the absolute angles 9^ of the segments. 
In order to be meaningful, such a function must be pe- 
riodic in each of those angles, Q{...,9^ + 27r, . . .) — 
Q{. . . ,9/^1, . . .), which in turn implies that it may be rep- 
resented by a discrete Fourier series: 



1, . . . ,Pjv) 



En i{p,iei+...+fj.i, 
^ ^ii,...,P,n'" 



(12) 



where integration yields the Fourier coefficient Qpi,...,;^^, 



2tt 



2tt 



X e-'(^i^~i+-+'^"^")p(6ii,...,0Ar) . (13) 

(It proves convenient here to employ a Fourier expansion 
scheme that differs from the angular decomposition (0) in 
its sign convention and normalization.) By combination 
of (|ll|l and (|13|l . we find the average value for the func- 
tional {Q{9i, . . . , 6'jv)) in terms of Fourier coefficients: 



)j ~ ^ ] Q tJ.l,...,flNPfiN ■ ■ • • ■ P/JJV + ... + A«l 



-V2,V2—l'5,---,l'N-l—VN,l'N 



^ PvnPvn-1 ■ • ■ ■ ■ Pvl 



(14) 
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In the second line of this equation, we employed an al- 
ternative indexing scheme {vi, . . . , j/jv} defined via: 

l^k = ^-lk + ■ ■ ■ + fJ-N < — > Mfc = t^A; - i^k+i , (15) 

(fiN = i^n)- Its advantages will become obvious below. 

We are particularly interested in the special case 
where the functional Q{9i, . . . ,9^) reduces to a prod- 
uct of N equivalent factors q{9i,): Q{6i,...,9n) = 
Yivli^i')- Under these circumstances, the Fourier com- 
ponents will factorize as well, and we have 
= <l^n ■ ■■■■Q^iN, where: 

q{9) = V e'^'q^ ^ q^ ^ ^ (^9 e-'^\{e) . 

(16) 

We now exploit these properties in order to simplify the 
average (|14|l . Upon inserting , we rearrange the sum- 
mations involved: 

UM - 1 

X . . . • ^ ^ qi,2-V3Pv2 ^ ^ 1ui-V2Pvi ■ (1''') 
1/2 fl 

Thus, we may rewrite {Q{9i, . . . , On)) in the form of a 
trace over the TV.th power of a matrix Z: 

(q(0i, . . . , 0n)) = E [^''lo.. = ^0 • 2^ . n , (18) 
1^1 

where [eo]k = Sok denotes a unit vector, and n = 
(...,1,1,1,.. .)-^ is composed of unit entries. The ma- 
trix Z itself is defined via: 

Zjk = qk-jPk , (19) 

and may be interpreted as the product of a circulant ma- 
trix Qjk — qk-j with a diagonal matrix Vjk = Pj^jk- 
(Occasionally, it is helpful to split off a scaling factor 
A from the coefficients q^: q^ = A'^g^, e. g. to ac- 
comodate functionals in a rotated frame of reference 
Qa{9i, . . . , 9n) = Q{9i - a, . . . ,0n - a), where = 
e"^"gi/ (|16|l . A brief inspection of p7|) shows that a slight 
generalization of p8|l will cover this extension: 

(Qa(0~i,...,0~^))=^A''^[Z^]„^^ .) (20) 

The expressions (|18|l and (|2U|I become especially handy 
when only a finite number of angular Fourier components 
Pk in @ are non-vanishing. If = for all \v\ > to, 
only the coefficients q-2m, ■ • ■ , q2m will enter the calcula- 
tion (irrespective of the nature of q{9)), and the problem 
in l(TH|) is effectively reduced to finding the iV.th power 
of a (2m -I- 1) X (2to -I- l)-matrix. Eigenvalue methods 
then offer an efficient algorithm to calculate the average 
for any N. In conclusion, we remark that our matrix for- 
mulation of the planar persistent random walk problem 



closely resembles the treatment of polymer chains in the 
rotational isomeric state (RlS)-approximation jTll |. and 
bears similarity to the Ising model of a one-dimensional 
ferromagnetic system The mathematical aspects of 
the representation of averages in Markoff chains through 
powers of a matrix have been summarized by MontroU 
'2(?| who incidentally applied the method to discuss the 
probability distribution function for a persistent random 
walk in a single dimension. 

As an important example for this formalism we exam- 
ine the characteristic function for the end-to-end prob- 
ability distribution c{uj,x) = (exp(ia; • 'L))n (where 
uj = (a;cosx,ti;sinx)"^). The projection of the vector 
L = {Lx,Ly)'^ onto the axis u) is given by a; • L = 
ujI J2k cos(6'fe — x)- Therefore, the function q{9) Hl()|) reads 
in this case: 

q{9) = e''^''=°''(^~-x) = ^i'^e-"'^J^(cjZ)e'''^^ , (21) 

V 

since q(9) coincides with the generating function for the 
Bessel functions of integer index 0| . We take advantage 
of the scaling symmetry (|20|l to identify A = ie~''^ and 
qv — Jv{^l), and consequently obtain for the character- 
istic function: 

c(c.,x) = Ei"'e"''''[^M)^]o. - (22) 

where Z(ujl)jk = Jk-j{y!l)Pk Ijltjl) . (In form of a recur- 
sive relation, this result is also stated in Ref. 0.) Equa- 
tion 1221) may be interpreted as the angular Fourier series 
for the characteristic function c(a;,x). Remarkably, if Pm 
is the non- vanishing Fourier component of largest index 
|m| in the series forp(0) then also the angular Fourier 
expansion of c{uj^x) will be a finite series limited to the 
range \v\ < \m\. In this sense, the angular decompo- 
sition of the end-to-end distribution directly reflects the 
properties of the angular bias p{9) between the individual 
segments that make up the walk. 

Finally, we note that our treatment covers Pearson's 
problem Q as the isotropic case p{9) = l/27r, or pk = 
Sok Q- Then, the characteristic function H22() is also 
isotropic and reduces to a simple power of a Bessel func- 
tion, cp(w) = Muil)^ 0. 

C. The probability distribution function 

It has been pointed out by various authors 0, that 
the probability distribution function w{R, (f>) = ((5[R — 
L(0i, . . . , 9n)])n is connected to the characteristic fimc- 
tion c(w, x) = (exp(ia; • J^))n via a Fourier transform in 
a;-space: 

w{R, <!>) = ^ J d'u: e---^ c(c., x) . (23) 

Since a; • R = ti;i?cos(x — 4>)i we may expand the ex- 
ponential in H23|) into the series of Bessel functions l|21|l . 
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Inserting the matrix- type expression for ciuj, x) ren- 
ders the integration over the polar angle x in trivial, 
and one obtains an angular Fourier series for the end-to- 
end distribution function: 

We may interpret this result as the generalization of 
Kluyver's formula |l7| towards the persistent random 
walk, to which it reduces in the isotropic case pk = Sok- 

1 f°° 

wp{R)^— ujdujJo{ujR)Jo{ujl)^ . (25) 
Jo 

As for c(u;,x), we again infer that the range of non- 
vanishing angular Fourier components in w{R, (p) equals 
the corresponding range for the local weight function p{9) 
concerning the individual angles between steps. 

While a rather elegant expression, H24|l is nevertheless 
not well suited towards the accurate determination of 
w(i?, (/)), since this pursuit involves the numerical eval- 
uation of an oscillating integrand (a sum over products 
of + 1 individual Bessel functions, see H17|l 'l. Indeed, 
since u;(ii!, 0) = for R> Nl [a. fact that we will take ad- 
vantage below), the distribution is not an analytic func- 
tion of R, and the asymptotics of the integrand in H24() 
suggests that w{R, (j)) shows singular behavior whenever 
the lengths ±i?, ±1 involved sum up to zero [ij. (For 
identical step lengths /, this occurs at regular intervals 
R = {N -2k)l, where k = 0, 1, 2, . . . , iV.) From the 
practical point of view, the eigenfunction expansion of 
w{R, (p) that we are going to develop now is much su- 
perior. On the other hand, the analytical expressions 
for the integral (|25() known for A'' < 3 provide a useful 
testing ground for any numerical attempt to evaluate the 
probability distribution function ( Section llllj) . 



D. The eigenfunction method 

Since the end-to-end distance i? of a walk of A' steps, 
regardless of its details, cannot exceed the contour length 
Nl of the chain, the distribution function w{R, (p) van- 
ishes outside the disc R < Nl. This renders a technique 
attractive that is familiar from the theory of field equa- 
tions under boundary conditions, viz., the expansion into 
eigenf unctions of a self-adjoint operator. For a function 
confined within a circle, the modes of a drum, i. c. solu- 
tions of the equation (A -I- A^^,)/mfc(^, 0) = 0, provide 
a natural choice. Properly normalized, these functions 
read HI: 



frnk{R,4') 



'.-'""^Jr^riZmkR/Nl) 

Nly/n Jm+i{zmk) 



(26) 



where Jm{z) is the Bessel function of order m, m G Z, 
Zmk denotes the k.th zero of Jm{z) (fc > 1), and Xmk = 
Zmk/Nl is the assigned eigenvalue. The functions H2t)|) 



form an orthonormal set that is furthermore complete, 
i. e., with = Lcosip we find for all R,L < Nl: 



oo oo 



m— — oo k—1 



Now, the probability distribution function 'w{R, (p) = 
((5[R— L(0i, . . . , 9n)\)n presents just the average of JSU, 
and we thus formally obtain 16]: 



oo oo 



W{R, (p) 



\ - sr^ Fjnk^ ^™'^Jm{ZmkR/Nl) 

^ 2^ ^ATi\2- T ^2 (28j 



as the Fourier-Bessel expansion of w{R,<p), where the 
weight coefficients Fmk are given by the averages: 



mk — \ *^ ^77 



N 



(29) 

The actual calculation of the coefficients Fmk fol- 
lows the general line of argument presented in Sec- 
tion III Bl where only a slight modification is required. 
The weights F^k present averages of functionals of the 
absolute angles of the segments 6i, in the sense of (|12|l . 
and we thus seek to establish the angular Fourier se- 
ries of Qmki9i, ...,eN) = e'^^Jm (zmkL/Nl). Although 
this may appear as a formidable task at first, an easily 
proven "projection theorem" ljAl|l for Bessel functions 
(Appendix^ will provide the desired result: 

QlT'Uiv = Sm.^,,+...+f,NJf,,{zmk/N) • . . . • Jf^^izmk/N) . 

(30) 

Inserting this result into 1141) and switching to the in- 
dexing scheme {vi, . . . , uj^} IjlSI) . we find that vi = to. 
Rearranging the remaining summations in the spirit of 
((T7|l yields at once: 

> J / Zmk \ J 



I/JV 



VI \ ^ ) . (31) 



N 



By comparison with JTSJ), we infer that this multiple sum 
can be recast again into the A^.th power of a matrix 
Z{zmk/N): Fmk = [Z{zmk/N)'^]^^, whose matrix el- 
ements are given by: 

[Z{zmk/N)]^^ = J^^,[^)p,. (32) 

In conjunction with (|28|l . this finally yields the Fourier- 
Bessel expansion of the end-to-end distribution function 
w{R, <p) of the persistent random walk: 



oo oo 



^oo^l (^0^7rJ„+i(z™fe)2 



X e 



-\m4> T I ^mkR 



(33) 
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Several remarks are in order here. First, we point out 
that the expansion coefficient F,nj. vanishes if the 
corresponding internal angular Fourier component pm 
does. Thus, the angular character of the end-to-end dis- 
tribution function w{R, (p) and the angular bias p{0) at 
the joints of individual segments (O are directly related. 
Like in Section III Bl a terminating Fourier expansion of 
p{9) implies a finite-size matrix Z{zmk/ N) in the deter- 
mination of the coefRcients in (|33|l . In particular, for 
Pearson's isotropic walk the radial distribution function 
wp{R) simphfies to: 



(34) 



a result first obtained by Barakat |3| . 

Second, series similar in structure to H33|l account for 
arbitrary derivatives of the probability distribution func- 
tion w(R, (p). We invoke the following differentiation re- 
lation 21] for the eigenfunctions fmk{R,(f>) I26|l : 



— ±i — 

dx dy 



V 'WJ 



"I 



-i(m^n) 



Nl 



'J, 



m=pn 



Nl 



(35) 



(which immediately yields the eigenfunction property 
Afrak{R,4') = -Knkfmk{R,(p) of thcsc solutions to the 
wave equation) and apply it to H33(l : 



(5, + idyY (d, - idy f wiR, <j>) = (-1)" 



oo 



m — — oo k—1 



X g-i(m+«-i-)</' J 



{Nlf^J^+l{ZrakY 



m-f-n — r 



Nl 



(36) 



All derivatives of the form df^dyW{R, (f) are then available 
as properly chosen linear combinations of the Fourier- 
Bessel series ()36|) . where n+r — The representation 

(|36|l is particularly useful in the random walk model of 
polymers .11], where the gradient V log w(i?, 0) is related 
to the average force exerted by the polymer chain. 

Finally, also all physically sensible moments {x^y'^)N 
[pL^v £ No) of the probability distribution function in 
the persistent walk can be represented as fairly simple 
sums involving the zeroes Zmk of Bessel functions ()B5|I . 
(|B6|I . Since the evaluation of these moments involves an 
integral over Bessel functions that is not contained in 
standard tables [3 1 the details of the calculation have 
been assembled in Appendix IbI 



III. EXAMPLE CALCULATIONS 

In the following, we probe different practical aspects 
of the Fourier-Bessel expansion H33|l derived in the pre- 
vious section. The numerical convergence of the method 



is the first topic to be examined in detail, and we present 
analytical and numerical results obtained on Pearson's 
original walk problem below (Section |lII Ap . Some of the 
distinctive properties of persistent walks arc highlighted 
in Section IlII Bl where we will study how the underlying 
angular bias p{9) infiuences the shape of the probability 
distribution function. Our examples include a "chiral" 
walk problem ^] that is biased towards left or right turns, 
and a model illustrating the transition from diffusive to 
elastic behavior as the stiffness of the chain is increased. 



A. The Pearson random walk 

Obviously, particularly simple results are obtained 
when the Fourier-Bessel method of Section IIIDI is ap- 
plied towards the isotropic Pearson random walk with 
p{0) ~ 1/2tt or pu = Sqi, 10). The angular sum in H33() is 
then absent, and the matrix power in the Fourier coeffi- 
cient reduces to a simple power, which allows convenient 
access to error estimates for the truncated sums necessar- 
ily encountered in a numerical evaluation of wp{R) via 
the expansion (|34l) . 

Unfortunately, the probability distribution function 
wp{R) for a random walk in two dimensions, unlike its 
three-dimensional counterpart 111 |2^ . generally defies 
evaluation in closed form. However, analytical expres- 
sions can be found for the shortest chains with N ^ 2 
and = 3. A formula for the two-segment chain distri- 
bution is available for arbitrary persistent walks (see Sec- 
tion CiTbIi, and reads here wp{R) = [ir^ Ry/AP - R'^]'^ 
(|45|l . For N = 3, the integral representation H25(l can 
be evaluated in closed form 2lj. Denoting the complete 
elliptic integral of the first kind by K{k), we find: 



wpiR) 



TT^{1 + R)y/{3l-R){l + R) 
Al I 



X K 



IR 



l + RV {il-R){l + R) 



(37) 



in the interval < R < L and 



wp{R) 



1 



2t:H^/IR 



K 



Al 



l + R {3l-R){l + R) 



IR 



(38) 

for / < i? < 3/. Note that wp{R) possesses a logarithmic 
singularity at i? = / and a steplike discontinuity at the 
maximum extension R — 31, where it approaches the 
value wp{3l) = 1 / {AVStt^ l'^) . (The division of wp{R) 
into distinct intervals of length 21 was already mentioned 
at the end of Section FlI CI ) Its exact shape is displayed 
in Figure 121 

These formulas provide a useful testing ground for the 
Fourier-Bessel expansion (|34|l . Clearly, this series can- 
not be convergent in a strict sense since wp{R) diverges 
as i? — > L For an approximate analysis, we estimate 
the asymptotic value of the sum terms Sp,{N,R) with 
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End-to-end distance R 

FIG. 2: The probability distribution function ■wp{R) for a 
three-step Pearson walk. The exact solution 1371 - 13811 (bold 
line) is shown together with the Fourier-Bessel expansion 
1341 1. truncated after fc = 10 terms (dashed line) and k = 100 
terms (solid line). For k > 1000, the approximation becomes 
visually indistinguishable from the exact solution. 




End-to-end distance R 

FIG. 3: Relative error erc\{k,R) incurred in the truncation 
of the Fourier-Bessel series 1341 for A'' = 3. From top to 
bottom, the solid lines display the relative error of the series 
for k — 100, k — 10*, and k = 10^, as compared to the exact 
result 1371 and 1381 . The heuristic error estimate 1401 (dashed 
lines) shows that the error scales inversely with k. 



large index /i. Since zq/j ~ '^{f^ ~ i) ^'^d Jo(Azo/t) ~ 
•\/2/7rAzo^ cos(Azo^ — 7r/4) holds as /i ^ oo [2l|, we find 
for i?> M/m: 



1 



7r(4^ - N 



2N 

-1) 



4iV 




(39) 



Thus, the series terms overall decay like Sfj,{N,R) ^ 
^-{N-i)/2 ^ Equation H39|l allows us to examine the ab- 
solute error e(fc, R) introduced by truncating the series 
(|34|l to terms with index /i < fc. This quantity is the sum 
of the neglected terms, e(fc,i?) — J2ii>k ^i^i-^^-^)' ^^'^ 
its asymptotics depend on the behavior of the oscillating 
terms in the second line of (|39|l . li N~ R/l is not en even 
integer, the average of the oscillating term vanishes, the 
series is of the alternating type, and for an estimate of 
e(fc, R) we may multiply Sk{N, R) with the largest period 
present in the product of cosine functions. These periods 
turn out to be proportional to the inverse distance of R 
to the interval boundaries {N — 2v)l (see Section Fll Cll . 
Indeed, for = 3 the maximum error e(fc, R) incurred in 
the truncation of the series is approximated well by the 
empiric formula: 



e{k,R) 



0.0076 /I 
klVlR V 2 



3/ 



21 



\l-R\ 31 -R 



(40) 



In the case that Nl—R is an even multiple of /, the Fourier 
series contains a strictly positive (non-oscillating) part 
that dominates the remainder. Its summation yields the 
estimate e(fc, R) ~ fc"^^"'^)/^ for the interval boundaries. 



except for iV = 3, where it reproduces the logarithmic 
singularity at i? = ^ H37|) . Examples for the numerical 
convergence of the Bessel-Fourier series for iV = 3 are 
shown in Figure Inland Figure |31 where the relative error 
erci(fc, R) = e(fc, R)/wp{R) is displayed together with the 
estimate (|40|) for various values of fc. 

Clearly, for fixed N the absolute convergence of the 
expansion (|34|l is slowest at the origin R = 0. We may 
then use arguments similar to H39|) to establish an abso- 
lute error bound e(fc) for the Fourier-Bessel series, which 
then yields fc as a function of the desired accuracy e: 



fc(6) 



AN 



i6r( 



4 ) 



7r4(iV-4)r(^] 



(41) 



Note that for short chains, the absolute series conver- 
gence at i? = (fc £-2/(Af-4)-j much worse than 



its typical convergence (fc ^ e 



-2/(N-l) 



). In particular. 



H41|l predicts the divergence of wp{R) for iV = 4 as 
R approaches the origin. Invoking the general relation 
wp{N + l,i? = 0) = wp{N,R = I) that follows at once 
from (|25|l . the analytic expressions H37() . (|38|1 immedi- 
ately confirm this assertion. With increasing number of 
segments TV, the number of sum terms fc(e) drops rapidly, 
and finally becomes of the order of the chain length it- 
self, in which case the estimate (|41|) becomes too conser- 
vative. In the long-chain regime 7rfc(e) ^ N the argu- 
ment of the Bessel function in 1341) tends towards zero, 



t2\N 



and the approximation Jo(zofe/^) ~ (l~-2^ofc/4-^ ) 
exp(— ZQJ./47V) instantly yields the asymptotic behavior: 



fc(e) ^y-Nlog{T:Ne) 



(42) 
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FIG. 4: Number of terms k in the Fourier-Bessel sum 13411 
required to guarantee a prescribed absolute accuracy e in 
wp{R). Solid lines, from top to bottom: e — 10~^^, 10~^, 
10~^. The dotted lines display the short-chain and long-chain 
approximations Mill . (1421 for e = 10"^'^, respectively. (Note 
that for iV > 316, wp{R) < 10"^ holds. For e = 10"^^ and 
N = 5,6, numerical convergence has not been achieved.) 

Therefore, for long chains, the numerical complexity of 
the Fourier-Bessel series approach grows merely with 
{N logN)^/^ . Figure 01 depicts numerical results for k(e) 
for various values of the absolute accuracy e, together 
with the estimates (|^ and . We conclude that these 
expressions closely reproduce fc(e) within their respective 
regimes of validity. 

Following this excursion on the numerical aspects of 
the Fourier-Bessel approach, we next study the result- 
ing probability distributions wp{R) (|34|) for the Pearson 
random walk. As a random process, wp{R) must be nor- 
mally distributed in the limit N ^ oo, and thus is com- 
pletely determined by the diffusion constant Dp — P 
(|10|l . After proper normalization, wp{R) asymptotically 
behaves like: 

-Hi^)-^exp(-^) , (43) 

assuming that R <^ Nl. For short walks of less than 
= 10 segments, deviations from the Gaussian limit be- 
come clearly apparent, and the separation of 'wp{R) into 
distinct functions for every interval of length 21 reveals 
itself. The nonanalytic behavior ofwp{R) at the interval 
boundaries R = (A^ — 2i/)Z is manifested by the presence 
of divergences (A^ = 3,4) and kinks (A^ = 4,5,6) in the 
distribution functions, as shown in Figure |S| 

Finally, we examine the force function Fn{R) as- 
signed to the random walk problem. Persistent ran- 
dom walks have been proposed as a model for the 
configuration of polymer chains this statis- 

tical interpretation, the bias function p{6) |(2J) is re- 
lated to the angular potential between adjacent seg- 
ments U{0) via (properly normalized) Boltzmann fac- 




End-to-end distance R/NI 

FIG. 5: Dependence of the probability distribution function 
wp{R) on the number of steps A'^. From top to bottom, the 
distributions for A = 4, 5, 6, 8, and = 12 are displayed. 
The nonanalytic behavior of wp{R) at the interval boundaries 
R = [N — 2v)l is clearly distinguishable for short chains. (For 
N = 4,, wp{R) diverges logarithmically as i? 14111 . ) As 
increases, the diffusive character prevails. The inset shows 
the distribution for A = 20 (solid line), together with the 
Gaussian approximation 14311 (dotted line). 

tors p{9) oc exp[— [/(0)/fcBT] which renders the Pear- 
son walk equivalent to the uncorrelated freely joined 
chain (FJC) model of the polymer. For a fixed end po- 
sition R of the polymer chain, the molecule exerts an 
entropic, temperature-dependent force whose average is 
given by the logarithmic derivative of the distribution 
function (F(R)) = fci3TVlogw(R), and its ab- 

solute value reduces to Fn[R) = —kBTwp{R)' /wp{R) 
for the isotropic walk. Like the distribution function it- 
self, its derivatives are conveniently calculated using the 
Fourier-Bessel formalism (|36|l , and the results of this ef- 
fort are displayed in Figure for the two-dimensional 
FJC model. 

For short chains, the singularities of wp{R) at the in- 
terval boundaries R = {N — 2v)l confer a conspicious 
irregular shape to the corresponding force curves Fisi{R). 
As A^ increases, Fm{R) becomes smoother, and as a func- 
tion of the scaled end-to-end distance R/Nl, approaches 
a well-defined limit Foo{R)- Its analytic form may be 
obtained by considering the FJC in the Gibbs ensem- 
ble, where the force F is kept constant, and R is al- 
lowed to vary llj. Its partition function is given by 
Z{F) = (exp(— L • F/fcsT)), and the average position 
of the chain end (R) under the applied force F follows 
at once from (R(F))Ar = -fc^ T V l ogZ(F). Using the 
technique presented in Section III Bl the function Z{F) 
can be cast as the trace over a matrix power (|18|l for an 
arbitrary persistent walk. Indeed, Z(F) may be inter- 
preted as the characteristic function c{ijj) for imaginary 
argument a; — iF/ZcsT (|22|l . and thus reduces for the 
isotropic walk to Zp{F) = hUF/k nT) ^ , where Io{z) 
denotes the modified Bessel function ^ll. Hence, for the 
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FIG. 6: Average force Fn{R) exerted by a freely rotating 
chain as a function of the end-to-end distance R. The plot 
shows the logarithmic derivative of 'wp{R) for various chain 
lengths A^. In the diffusive limit N ^ oo, the force curves 
(solid lines, bottom to top: N = 8, 12, 20) tend towards the 
function Foo{R) ll^ (dotted line). The inset displays Fn{R) 
for short chains (solid line: A = 4, dashed line: N — 5, dotted 
line: A' = 6). The peculiar appearance of these functions is 
linked to the nonanalytic behavior of 'Wp{R) (Figure 0. 



FJC model we find (i?(F))jv = NlIi{u)/Io{u), where 
u — IF/ksT. In the limit — s- oo, we may identify 
the average position {R{F))n with the actual position 
R, and likewise Fn{R) with F. This allows to revert the 
relation to express Foo{R) as a function of R/Nl: 



Foo{R) 



knT 



U{R/Nl) , 



(44) 



where U{z) is the inverse function to Ii{z)/Iq{z). It is 
defined for jzj < 1 and behaves like U{z) — 2z + z'^ + 



5 ^5 



for small \z\, while its Laurent series near z = 1 



is given hy U{1 — z) = j-. 



2z 



B. Properties of persistent walks 

Analytical methods towards the calculation of the 
probability distribution function w{R, (p) for the more 
general persistent walk problem are usually limited to the 
combination of only two segments (A^ = 2). In this case, 
for a fixed endpoint L = (i? cos </>, R sin </>) only two possi- 
ble locations for the intermediate joint exist, and simple 
geometrical arguments lead to the following result: 



wiR, 0) 



i?V4/2 - i?2 

X [p(27)p(</>- 



7)+p(-27)p(0 + 7)] , (45) 



where we used 7 = arccos(i?/2Z) to denote the angle be- 
tween the end-to-end vector and the individual segments. 
(For the uncorrelated Pearson walk, the term in brackets 



FIG. 7: The elastic model of a stiff persistent random walk. 
The chains are replaced by circular arcs of fixed contour length 
Nl. When the curvature of the arc is changed, its end point 
moves on the geometrical curve R = Nl sin ( 



reduces to l/(27r^).) Likewise, the limit of very extended 
walks (A^ — > cxd) is amenable to an analytical treatment. 
Here, the end-to-end probability function w{R, (p) must 
be again normally distributed, and is completely deter- 



mined by the persistence lengths {L^) 
and the diffusion constant Dao l(Tn|) : 



u;(R) 



t:NDo 



■ exp 



(R- 



©, {Ly)oo ®, 



(46) 



We point out that the distribution w{R, cf)) for sufficiently 
large A^ depends merely on the single angular Fourier 
component pi . Little can be said a priori about the be- 
havior of w{R, (j)) for intermediate chain lengths A^. For- 
tunately, the Fourier-Bessel expansion provides an 
efficient means to determine the probability distribution 
function numerically. 

In the following, we select normalized angular bias 
functions p{0) of the general form: 



' TTV{2i + 1) 



(47) 



where ^ > and (3 are real parameters. The correspond- 
ing angular Fourier components p^, jSJ are given by |l8l | : 



r(e + l)2e'''^ 



(48) 



In practice, the coefficients Pu{£,,f3) are best calculated 
by the simple recursion formula pi^+i = {£, — i^)Pu / {S, + 
1^+1) {h' > 0), where normalization requires po = 1. We 
find in particular pi = e'''^/ (^ -f 1) and thus: 



l-2(e + l)sin2| 



D„ 



1 + 4C(e + 1) sin2 I V (e + l)sin/3 



(49) 
(50) 
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in (|46|l . Note that the angle /3 in (|47|l controls the chiral 
bias of the walk Q ; a symmetrical distribution will ensue 
if /3 = is chosen. (Then, l|49|l simplifies to {Lx)oo — l^, 
{Ly)oo — 0, and Doo = ^^(2^ + !)•) Finally, we remark 
that the series in H47|l terminates whenever ^ is an integer. 
(Otherwise, p{6) is not analytic at its zero 9o = tt + (3.) 

We start out with a study of the dependence of w{R, cf)) 
on the number of steps N for the fixed angular weight 
p{9) = 3I cos4[i(6' - f )], which is of the form gTj) with 
£, — 2 and /3 — 30°. The angular bias evidently favors 
a "spiral" configuration of the chain in which adjacent 
segments are joined at an angle of /3, which translates to 
an absolute angle = {N+l)P/2 of the end-to-end vector 
with respect to the x-axis. However, configurations that 
end up near the origin constitute the vast majority in 
the phase space {9i, . . . , On} of the walk. For sufficiently 
large N, the latter contribution always dominates and 
ultimately yields the Gaussian distribution H46(l . while 
the shape of wlR, (p) is governed by the intrinsic stiffness 
of the chain for short walks. 

In this regime, large deviations from the ideal interseg- 
ment angle [3 are suppressed by the angular bias function 
p{0)^ so walks that significantly contribute to w{R,(j)) 
show small variations of their bond angles 6k around (3. 
While these variations may significantly alter the total 
turning angle of the persistent walk, they usually have 
little effect on the end-to-end distance R, at least when 
the preferred configuration is mostly stretched (i. e., Nf] 
is small). The situation then resembles the elastic re- 
sponse of a thin rod (or, in the case of chiral bias /3 7^ 0, 
a circular spring), with one end fixed at the origin: While 
easily yielding to bending forces, the rod (spring) is re- 
silient to compression, and resists changes of its length. 
Under applied force, the condition of constant curva- 
ture thus requires that the rod bends into a circular arc 
of fixed length A, and a simple geometrical argument 
shows that the free end of the rod moves on the curve 
R = A sin 0/0. The shape of this curve is displayed in 
Figure [7| 

A numerical calculation of the probability distribution 
function w{R, (p) for the aforementioned angular hias p{9) 
confirms these deliberations. In Figure |S1 we present 
w{R, 0) for walks of length = 4, iV = 7, and TV = 10. 
For the short walk, 'w{R^ 0) is concentrated along a nar- 
row groove whose shape is reminiscent of the geometri- 
cal curve presented in Figure while for the extended 
chain, the distribution is largely of Gaussian shape, al- 
beit shifted from the origin, as predicted in H49|) . In the 
transitory case N = 7, the end-to-end distribution fimc- 
tion conspicuously combines the typical features of both 
limiting cases, as both the central hump and a protruding 
ridge are clearly visible. (Note that the isolated spike in 
the center of the left image is not a numerical artifact; the 
actual distribution for A^ = 4, like the corresponding un- 
correlated Pearson walk (Figure^ contains a logarithmic 
singularity at the origin. In the calculation on display, 
we chose to truncate the Fourier-Bessel series (I33|l after 
k = 10'* terms.) 



In a second simulation, we investigate the probabil- 
ity distribution function w{R, (f) in a persistent walk 
of A^ = 20 segments as a function of the stiffness of 
the joints, i. e., the parameter ^ in the angular bias 
function p{9) (|47l) . In this example, we chose a sym- 
metric distribution (/3 = 0) and calculated the aver- 
age entropic force (F(R)) = fcsT| V logz«(R)| (see Sec- 
tion nil A|l for the three stiffness parameters ^ = 2, 
^ = 4, and ^ = 6. (Under these circumstances, con- 
vergence of the Fourier-Bessel series H33f) and (|36|) was 
achieved after summing k — 250 terms.) Lines of con- 
stant force are plotted in Figure El For the softest chain 
(^ — 2), the pattern resembles a set of equidistant, 
concentric rings that has been shifted from the origin. 
This observation implies that the corresponding distri- 
bution z«(_R, (f) has approached the long-chain limit 14t)|) 
characterized by an isotropic harmonic force field i^(R) 
directed towards a minimum Rmin of the free energy 
t/(R) — — fcsnogw(R). Its position is given by the per- 
sistence length vector Rmin = (L)oo = (Z^,0)^ (ggi, and 
the elastic constants Cjk = 9j9fcC/(Rmin) there read Cxx — 
Cyy = 2kBT/ND^ = 2kBT/[{2i + l)Nl^], and e,y = 
0. The numerically determined values Xmin = 2.3224/, 
e^^ = 1.522 • 10-2 keT/P, and eyy = 1.438 • lO'^ keT/P 
still deviate considerably from their counterparts in the 
diffusive limit, a^min = 21 and exx = ^yy = 0.02 ksT/P , 
however. 

The force distributions obtained for the stiffer walks 
with ^ = 4 and ^ = 6 bear little resemblance with 
the diffusive circular pattern observed for ^ = 2. For 
^ = 4 (center panel in Figure |5J), the force field is dis- 
tinctively flat over a wide range of end-to-end vectors 
R, and the free energy U{IV) actually features two local 
minima x< = 3.9808/, x> = 15.160/ on the x-axis. The 
latter represents the maximum of the probability den- 
sity with w{x^,0) « 1.91 • 10~^l~^, but their difference 
in free energy is only AU — 0.2336 fesT. The minima 
are separated by a shallow maximum near Xmax = 9.15/ 
{AU = 0.28A6 ksT). It is instructive to compare the 
elastic constants ejk at both minima. For x<, we find 
an almost isotropic response {exx — 7.421 • 10~^ ksT/P, 
Cyy = 7.009 • 10^^ ksT/P), while the free energy surface 
near a;> is strongly prolate {exx = 9.529 • 10~^ ksT/P, 
Cyy = 2.557 • 10^^ ksT/P). Clearly, this anisotropy indi- 
cates the onset of elastic behavior. 

For f = 6 (right panel), elastic response becomes the 
dominant feature of the distribution which is predomi- 
nantly concentrated along a narrow arc, while a rather 
fiat plateau in its center still persists. The maximum of 
w{R, (p) has shifted further out to a;> = 16.790 /. (An ex- 
tremely shallow local 2.4924/ {AU ^ 
1.222 fcsT) with a barrier of only 0.009 /csT is contained 
within the plateau structure.) The anisotropy of the free 
energy surface near x> has become even more marked, 
where the bending mode {eyy — 3.232 • 10"'^ ksT/P) is 
almost 75 times softer than the response towards com- 
pression {txx = 2.405 • IQ-^ ksT/P). Empirically, the 
shape of the trough structure in U (R) is accurately re- 
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FIG. 8: Probability density distribution w{R,<j)) for a persistent random walk with p{9) = ^cos'*[|(6 — 30°)], made up of 
= 4 (left), N — 7 (center), and A'' = 10 (right) segments. The transition from a stiff chain with a narrow distribution to a 
diffusive structure with Gaussian character is obvious. In the intermediate regime, elements of both limiting cases coexist. 
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FIG. 9: Lines of constant force F{R, 0) for persistent walks composed of A*" = 20 segments. The angular correlation is of the 
form p{d) oc cos^^(6'/2) 1471 1. With increasing forward bias (left to right: ^ = 2, 4, 6), a transition from largely diffusive behavior 
(a pattern of equidistant concentric circles) to the elastic response of a stiff rod (that is easily bent but resists compression) is 
observed. The dashed line in the right plot has been adapted from Figure |7| (Spacing of force lines: 0.03 ksT /I.) 



produced by the elastic spring model (with A = Nl) put 
forward in Figure [3 if we allow for a displacement of 
its origin so that the curve R{(l)) = Nl sin (p/ip passes 
through xy (dashed line in Figure EJ. The elastic com- 
pressive response is limited to small deviations from the 
equilibrium point a;>: For applied loads greater than a 
critical force i^crit = 0.195 ksT/l (at Xcrit = 13.90 Z), the 
stretched chain will buckle and collapse into the more 
globular shapes making up the plateau. This again re- 
sembles the behavior of an elastic rod under stress 1231. 



IV. DISCUSSION 

In this article, we examined persistent random walks in 
two dimensions, where the relative orientation of succes- 
sive segments is correlated through an angular bias func- 
tion p{0). We found that the eigenfunction expansion 
into a Fourier-Bessel series provides a convenient and 
efficient means to calculate the end-to-end probability 
distribution function w(R) encountered in a planar ran- 
dom walk except for extremely short chains. Related se- 
ries expansions also cover all derivatives and moments of 
u'(R). We used Pearson's original random walk problem 
to assess the numerical aspects of the Fourier-Bessel se- 
ries and found that the number of sum terms required to 



achieve a desired accuracy initially sharply decreases with 
the chain length N and, for very extended walks, grows 
like (TV log A^)^/^. This property renders the method 
attractive in studying the properties of chains of inter- 
mediate length whose distribution function i/7(R) may 
strongly deviate from the diffusive behavior expected in 
the limit N oo. Even for uncorrelated planar walks, 
the singularities contained in wp{R) cause conspicuous 
features in the probability distribution and force curves 
that only gradually fade away as A'' increases. Naturally, 
correlated walks generate a much richer set of phenomena 
that we only have begun to explore in this work. Here, 
we briefly studied "chiral" walks with asymmetric angu- 
lar bias p{9) that have a tendency to loop, and inquired 
into the dependence of w(R) on the stiffness of the in- 
tersegment joints. Clear evidence of a transition from 
diffusive character to a form of elastic response familiar 
from the elastomechanics of a stiff rod, was obtained. 

The latter observation raises questions that deserve 
further examination. It appears worthwhile to study the 
dependence of the crossover from diffusive to elastic be- 
havior as a function of N and the bias p{0). Character- 
istic parameters of the distribution, like the location of 
the minima Rmin of the free energy, the elastic moduli 
ejk and the critical force -Fcrit, may obey scaling laws. 
The existence of several minima of ?7(R) suggests that 
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the average extension of a stiff chain under an apphed ex- 
ternal force shows hysteresis. The formahsm developed 
in this article directly applies to these problems. 

Clearly, for many physical applications, e. g., in poly- 
mer theory [ij, realistic models will require the random 
walk to take place in three spatial dimensions. While the 
analysis of the uncorrelated random walk problem en- 
counters no major obstacles 15, 23], the extension of the 
eigenfunction expansion technique to persistent walks in 
three-dimensional space presents a challenging problem, 
since the angular bias p{9, (p) then generally will describe 
the relative orientation of a segment and its two pre- 
decessors, unless p{9, </>) is independent of the dihedral 
angle 4>. Alternative methods, like the transfer operator 
formalism , may prove to be better suited for the nu- 
merical determination of the probability density function 
w(R) in this case. 
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perform this projection also in an alternative reference 
frame which is rotated by the angle a, i. e. with respect 
to an arbitrary direction. Since this operation does not 
affect the lengths involved and shifts the angles (p and 
9^ merely by the constant value a, wc find the relation 
_Rcos((/) — a) = li, cos(0jy — a), or in exponential form: 



N 

i?cos(0— q) ^iulf^ cos{9u—a) 



i/=l 



(A2) 



for all u €E C and a S R. We employ the generat- 
ing function expansion for the Bessel functions J^{u), 
giucosa ^ i'=e"=" Jfc(u) (EU on both sides of with 
u = z/ A to obtain the equality: 



E(ie-)%-^j.(zi^). y: 



R 



le 



—ia\ ^i~\~---~\~^N 



Ai ,...,Aiv 

(A3) 

The Fourier decomposition of l|A3(l with respect to the 
angle a, i. e. an integration ^ da e''^" • • . ., then imme- 
diately yields the projection theorem (jAl|) . 



APPENDIX B: EVALUATION OF MOMENTS 



APPENDIX A: THE PROJECTION THEOREM 

This appendix provides a brief mathematical deriva- 
tion of the projection theorem for Bessel functions which 
forms the basis of the Fourier-Bessel expansions derived 
in Section lll Dl If {R, (/>) denotes the location of the chain 
end in polar coordinates for a walk consisting of N steps 
of individual lengths {h, . . . , In} characterized by abso- 
lute angles {^i, . . . ,9n} (see Figure Pi, this proposition 
states that: 



R 



E 



„i(Aiei+...+ANejv) 



Xl + ... + \N=fJ- 

X Jx,{zh/A) ■ Jx„ (zIn/A) , (Al) 

where A — li + . . . + 1^ represents the contour length 
of the walk, e Z is an integer, and z e C denotes an 
arbitrary complex number. (For the case of equal step 
lengths li, = I primarily considered in this article, the 
arguments of the Bessel functions simplify to zR/Nl on 
the l.h.s. and z/N throughout the r.h.s. of this equation, 
respectively.) 

We now sketch the proof of IjAip . From Figure ^ it 
is clear that the projection A^; of the end-to-end vec- 
tor on the X-axis is simply the sum of the correpond- 
ing projections of the steps that make up the walk: 
Rcoscj) = cos9„ = 1; • • ■ J It is convenient to 



The Fourier-Bessel series H33|l for ii;(i?, (f) permits fast 
and accurate numerical evaluation of all moments of the 
probability distribution in the persistent planar walk. 
For this two-dimensional problem, the moments are av- 
erages of the general form (L^L^')jv. Since it is easily 
verified that they may be expressed as a finite combi- 
nation of the complex averages (i?'"+^"e'™'^)Ar, where 
j + k — m + 2n and m,n £ Nq, we prefer to use the 
radial definition, as it leads to shorter expressions. Ex- 
amples are encountered in Section III Al where we calcu- 
lated the average position {L^ + iLy)^ = {Re^'^)N ifTTjl - 
corresponding to m = 1 and n = 0, and the mean square 



displacement (L^ -I- L'i 



>N 



{R'^)n (113, where to = 



and n — 1. (We also note in advance that for m = n = 0, 
conservation of probability demands {R'^)n = !■) 

The determination of the radial moments involves an 
integral over Bessel functions that apparently is not listed 
in the major tables ^T§\. Hence, we present a short sketch 
of its derivation. We start out with a recursion formula 
stated by Watson |2l|. Introducing the notation: 



(«) 



mk 



Au U 



m+2n+l 



Jra{u) , 



(Bl) 



where TO,n > are integers, and z^k denotes the fc.th 
zero of the Bessel function Jm{u), this recurrence reads: 



An) 



Z, 



m+2n+l 



Jm+i{zmk) - 4n(TO + n)4"fc . (B2) 
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Since the recursion terminates for n = 0, we may add 
up the terms and rearrange them as the product of 
Jm+i{zmk) with a polynomial in the zero Zmk'- 



An) _ J (, -.Sr (-4)" "n!(m + n)! +2.+1 



i/=0 



(B3) 

From the series expansion H33|l . we now calculate the 
radial moments ^/j™+2ngim0^ ^ ^ Since the angular inte- 
gration is trivial, using the abbreviation (|B1|) in the first 
step we obtain: 

{■■■)n= di?i?™+2"+M d(Pe'""t'w{R,(P) 
Jo Jo 

°° \Z(z , /N)^] /"'"^ 

m+2n+2 j / X2 ' 
/c=l nifc Jm+iyZmk) 

where the elements of the matrix Z are given by Zji — 
Ji-j{Zmk/N)pi (see Section III D|) . Using (|B3ll . we rear- 
range this sum into a series over powers of Zmk- It is 
convenient to define a set of auxiliary sums Bmr involv- 
ing the zeroes of the Bessel functions Jm{u): 



^ [Z{z^k/N)^] 



Om 



mr / , 2r+l J / \ 
k=l ^mk ■Jra+l\Zmk) 



(B5) 



The final result for the radial moments then can be cast 
in the following form: 



m-\-2n 



(-l)"-n!(m + n)! (^) 

v=0 ^ ' 



shows. In the calculation of Bmr (|B5ll . it is thus conve- 
nient to set = Q for > m + r.) 

Finally, we present several examples to illustrate our 
result for the moments of w{R,<j)). For n — 0, the se- 
ries (|B6ll reduces to a single term: {{L^ + iLy)"'^)j\[ = 
2{Niy"B^Q . In the case to = 0, conservation of prob- 
ability then requires -Bq^'' = 1/2, irrespective of N and 
p{0). Similarly, comparison with Q yields for m = 1: 



B 



(N) 
10 



Pl(l-Pf) 

2N{l-pi) 



(B7) 



We may use these relations to keep track of the accu- 
racy in the summation of the Fourier-Bessel series (|33ll . 
Similarly, from ljB6p we obtain all purely radial moments 



2n\ 



)N 



of even order. They read for n < 3: 





= {Nlf (1 
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(B8) 
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-32<)- 


hl28i?(^)) , (B9) 




^{Nlf (1 


-72<)- 


hll52i?W_ 4608sW) , 






(BIO) 



where we already inserted Bqq'' — 1/2. Clearly, an ex- 
pression for -Bq^'' in terms of the angular bias Fourier 
components p±i akin to (|B7p is available from our ear- 
lier result for the mean square displacement (|10|l . From 
there, we also infer that {R'^)n in the limit N 00 only 
grows linearly with N (|18l) . Hence, asymptotically the 
Bessel sum coefficient in (jB8|l approaches the value 
B^ = i. Similar considerations for the higher moments 
in conjunction with ljB5|) then show that: 



For not too small N, the series converge 
rapidly, thus permitting efficient numerical evaluation of 
the moments using (|B6|I . (For practical purposes, we note 
that the average (i?™+2»e'""^) ^ is a function of the 
angular Fourier coefhcients po,P±i, ■ ■ ■ ,P±(m.+n) only, as 
the method of direct evaluation exposed in Section FlI Al 



B, 



(00) 

Or 



00 

fe=l ^Ofc 



Mzok) 22'-+i(r!)2 



(Bll) 



where the coefficients a,, form a sequence of integers: 
flo = 1, fli = 1, 02 = 3, 03 = 19, etc. 
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